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ABSTRACT 

We performed a new analysis of B and R light curves of a sample of PG 
quasars. We confirm the variability-redshift correlation and its explanation in 
terms of spectral variability, coupled with the increase of rest-frame observing 
frequency for quasars at high redshift. The analysis of the instantaneous spectral 
slope for the whole quasar samples indicates both an inter-QSO and intra-QSO 
variability-luminosity correlation. Numerical simulations show that the latter 
correlation cannot be entirely due to the addition of the host galaxy emission 
to a nuclear spectrum of variable luminosity but constant shape, implying a 
spectral variability of the nuclear component. Changes of accretion rate are also 
insufficient to explain the amount of spectral variation, while hot spots possibly 
caused by local disk instabilities can explain the observations. 

Subject headings: galaxies: active - galaxies: photometry - galaxies: Seyfert - 
quasars: general - quasars : variability 



1. INTRODUCTION 



Although variability plays a key role in constraining the size of the central engine of 
active galactic nuclei, yet its physical origin remains substantially unknown. Even restricting 
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to the class of non-Blazar objects, the most diverse mechanism have been proposed in recent 
years, including gravitational lensing due to intervening matter (Hawkins 1996), supernovae 
explosions (Aretxaga, Cid Fcrnandes & Terlevich 1997), instabilities in the accretion disk 
(Kawaguchi et al. 1998) and star collisions (Torricelli-Ciamponi et al. 2000). For a small 
number of low redshift objects, a multi- wavelength monitoring with adequate time sampling 
and resolution allows the interpretation of changes of the spectral energy distribution (SED) 
in terms of an interplay of emission components with different spectral and variability prop- 
erties (see Ulrich, Maraschi, and Urry (1997) for a general review and Courvoisier (1998) 
for the case of 3C273). In the near infrared-optical- UV bands, variabihty studies indicate 
a hardening of the spectrum in the bright phase (Cutri et al. 1985; Edelson, Krolik & 
Pike 1990; Kinney et al. 1991; Paltani & Courvoisier 1994). So far, however, most of 
the statistical information on AGN variability, derives from single-band light curves of mag- 
nitude limited samples of objects (Angione & Smith 1972; Bonoli et al. 1979; Hawkins 
1983; Trcvcse et al. 1989; Cristiani, Vio & Andreani 1990; Trevese et al. 1994; Hook 
et al. 1994; Bershady, Trcvcse & Kron 1998). In the case of magnitude limited samples 
the analysis is complicated by the strong luminosity-redshift (L-z) correlation, caused by 
the crowding of objects towards the limiting flux. As a consequence it is difficult to isolate 
any intrinsic variability-luminosity (v-L) and variability-redshift (v-z) correlation. Moreover, 
the results of these analyses depend on the specific variability index adopted, as shown by 
Giallongo, Trevese & Vagnctti (1991) who found a positive v-z correlation through a vari- 
ability index defined on the basis of the rest-frame structure function. The existence of an 
average increase of variability with redshift was later confirmed by Cristiani et al. (1996). It 
turns out to be consistent with the suggestion of Giallongo, Trevese & Vagnetti (1991) that 
QSOs at high redshift appear more variable since they are observed at a higher rest-frame 
frequency, where the variability is stronger (hardening in the bright phase and vice versa). 
A direct statistical evidence of the spectral variability, in terms of an average change of the 
B-R color with the variation of the B magnitude, was found by Giveon et al. (1999) as 
part of a statistical analysis of the B and R light curves of a sample of PC QSOs. A similar 
evidence has been found by Trevese, Kron, & Bunone (2001) for the faint QSO sample of 
SA 57 (Koo, Kron & Cudworth 1986). Unfortunately, most (if not all) of the variability 
mechanisms proposed so far imply a hardening of optical-UV spectrum in the bright phase, 
so that they cannot be discriminated on purely qualitative grounds. As a first step towards 
constraining variability models quantitatively Trevese, Kron, & Bunone (2001) compared 
with the observations a simple model where both spectral slope changes and brightness vari- 
ations are due to temperature changes of an emitting black body. In the present paper, we 
present a new analysis of the B and R light curves of a sample of PG QSOs made available 
by Giveon et al. (1999), we analyze the consistency with previous results, and consider 
possible different sources of spectral variability. The paper is organized as follows: in section 
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2 we summarize the characteristics of the Giveon et al. (1999) data, in section 3 we analyze 
the v-z correlation, in section 4 we discuss the parameters adopted for the spectral variation 
analysis, in section 5 we consider the effect of the host-galaxy on SED variations, in section 6 
we compare the observed spectral variations with those produced by changes of the accretion 
rate, and discuss the consistency of the observations with a simple model consisting of hot 
spots on the accretion disk. Section 7 contains a summary and the conclusions. We adopt 
Ho — 50 km s~^ Mpc"-*^, qo = 0.5, unless otherwise stated. 



2. THE DATA 

The present analysis is based on the light curves made available to the community by 
the Wise Observatory group (Giveon ct al. 1999). The sample consist of 42 PG QSOs 
selected to be nearby, i.e. z < 0.4, and bright, i.e. B < 16 mag. 

The observations were made in the Johnson-Cousins B and R bands with the 1 m Wise 
Observatory telescope. The total duration of the campaign was 7 years and the median 
observing interval of the objects was 39 days. The r.m.s. photometric uncertainty is ~ 0.01 
mag and ~ 0.02 mag in the B and R bands respectively. We refer to the original paper of 
Giveon et al. (1999) for all details concerning observations, calibration etc. The full B and 
R light curves of the entire sample were made available in electronic form. 

Giveon et al. (1999) present the analysis of the correlation of different variability 
properties with other properties like luminosity, redshift, radio power, various line intensities 
and the X-ray spectral slope. In particular they show a correlation between the color changes 
A{B — R) and the brightness variations AB and AR, corresponding to an average hardening 
of the spectrum in the bright phase (and vice versa). They do not find a correlation between 
variability and redshift, at variance with Giallongo, Trevese & Vagnetti (1991), Trcvese et 
al. (1994), Cristiani et al. (1996), but they ascribe this to the difficulty of disentangling 
v-z, v-L and L-z correlations in the sample. We stress that this is made particularly difficult 
by the small redshift interval spanned by the sample. We perform a new analysis of both the 
v-z correlation and the spectral slope changes in the following sections. For these purposes 
we have dereddened the data, using the extinction calculator of the NASA Extragalactic 
Database (NED). 
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3. VARIABILITY-REDSHIFT CORRELATION 

To measure the amplitude of variability we define, for each object, the first order struc- 
ture function, as in Di Clemente et al. (1996) (but omitting the subscript "1"): 

S{t, At) = [{^\m{t + r)-m{t)f - <7^)]^ (1) 

where m{t) is either the B or the R magnitude, t is the rest-frame time, r is the time lag 
between the observations, cr„ is the relevant r.m.s. noise and the bar indicates the average 
taken over all the pairs of observations lying in the time interval r ± At. In this definition 
we adopt the square average of the absolute values of the difference instead of the average 
of the square difference as in Di Clemente et al. (1996), since the former quantity is less 
sensitive to outliers (Hook et al. 1994). The 7r/2 factor normalizes S to the r.m.s. value 
in the case of a Gaussian distribution. The adopted value of At is the result of a trade-off 
between time resolution and statistical uncertainty. 

In the following we define four variability indices Si{T ± At), with i — B,R, and 
T = 0.3 ± 0.09 yr, 2.0 ± 0.6 yr. The subscripts B and R refer to the observing band and the 
values of r and Ar have been chosen for comparison with previous analyses (Di Clemente 
et al. 1996). None of these four indices shows a significant correlation with redshift when 
the whole sample is considered, confirming the result of Giveon et al. (1999). 

However, to disentangle the v-L, v-z and L-z correlations, we can restrict the analysis 
to a magnitude bin —23.5 < Mb < —22.5, around the average absolute magnitude of the 
sample < Mb >= —22.75. The result is shown in Figure 1. In this case we find a v-z 
correlation coefficient r^^z = 0.39, which is marginally significant {P{> r) = 0.09) despite 
the small number of objects (19) in the bin. To examine the dependence of variability on 
redshift we take the ensemble averages of the four variability indices defined above, over the 
same subsample of 19 objects. For each observing band we compute the average rest-frame 
observing frequency of the sample. To compare the result with Di Clemente et al. (1996) 
we must take into account the dependence of variability on magnitude. For this purpose, 
we reduce S by an amount AS* = {dS/ OMb) AMb where for S{Mb) we adopt model A of 
Cristiani et al. (1996) and AMb is the difference between the average absolute magnitudes 
of the present sample and the sample of Di Clemente et al. (1996). Figure 2 shows the 
increase of variability with the observing frequency. At each frequency variability is larger 
for larger time lag r (due to the increase of the structure function from T = 0.3yrtOT = 2 
yr). The new points, obtained from the present analysis, are consistent with the previous 
results, which were obtained from different samples and observing frequencies. The general 
trend can be quantified as: 
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dSi/d\ogVrest = dSi/d\og{l + ~ 0.25 - 0.3 (2) 

and is consistent with Giallongo, Trevese & Vagnetti (1991), Cristiani et al. (1996) and 
Di Clemente et al. (1996). Therefore the increase of variabihty with redshift and its 
interpretation in terms of an average increase of the amphtude of variabihty with frequency 
are confirmed by the present results despite their poor statistical basis. 

Hawkins (1996) finds that the variability timcscalc is independent of redshift, i.e. it is 
not affected by cosmological time dilation. He suggests that the effect can be explained in 
the framework of gravitational microlensing caused by intervening matter. Hawkins (2001) 
estimates that a possible decrease of variability timescale with frequency is not sufficient 
to compensate for the cosmological time dilation. The frequency dependence of variability 
amplitude can also compensate, at least in part, for time dilation. We note that from Eq. 
2 we can estimate a vertical shift of ~ 0.1 between B and R power spectra, approximately 
consistent with the horizontal shift of ~ 0.06 found by Hawkins (2001) (see his Figure 3). 
However, a detailed evaluation of the amplitude shift of the variability power spectrum as a 
function of redshift would require the knowledge of the redshift, luminosity, and observing 
time distributions of the sample. 

We also stress that the frequency dependence of microlensing (Alexander 1995) should 
be compared with the results shown in Fig. 2 and Eq. 2. 

4. THE SPECTRAL SLOPE a AND THE SPECTRAL VARIABILITY 

PARAMETER p 

The average increase of variability with frequency can be interpreted in terms of a 
change of the spectral slope with luminosity. However the statistical evidence discussed 
in the previous section, and in Giallongo, Trevese & Vagnetti (1991), Di Clemente et al. 
(1996), is only indirect and must be confirmed by observations of individual objects in at 
least two bands. This was done by Giveon et al. (1999) in terms of correlation of color 
changes A{B — R) with brightness variations AB or AR and by Trevese, Kron, & Bunone 
(2001) in terms of changes of the spectral slope a, defined by L^, ~ i/", where Li, is the 
intrinsic power per unit frequency. The latter work is based on U,Bj,F, N observations 
at two epochs of the sample of QSOs of the Selected Area 57. The increase of spectral 
slope is computed as an average value over the ensemble of objects, and compared with a 
simple model consisting of a black body subject to small temperature changes. The data 
of Giveon et al. (1999), which contain on average 40 B and R observations of each QSO, 
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allow the statistical analysis of spectral slope changes of each object, for which we compute 
the instantaneous slope: 

a{t) ^ log(L.jL.J/log(.B/.«) = J-'^[iB-R)- {Bo-Ro)] _ 

where ttIo = —2.5 log fo are the zero points of B and R photometric bands respectively (Cox 
2000). Taking into account that typical variabihty time scales are Ri 1 yr (see e.g. Trevese 
et al. (1994)) we regard as simultaneous the observations within a time interval of 9 hours. 
In Figure 3 we report the instantaneous value of the spectral slope a as a function of the 
intrinsic luminosity L,^^. Each small cloud of points represents a single QSO in different 
luminosity states. The relevant regression lines are reported on each cloud. They show, 
with a few exceptions, a positive correlation between the spectral slope and the intrinsic 
luminosity, which we call intra-QSO a — L correlation. The distribution of clouds in Figure 
3 also shows that the average slope of each QSO tends to be larger for brighter objects, 
forming a sort of QSO main sequence in the a — L plane (Trevese & Vagnetti 2001a). This 
inter-QSO a — L correlation can be quantified considering for each QSO the average values 
(a) and (L^g) of the slope and luminosity. The correlation ra-L — 0.58 is highly significant: 
P(> r) = 6 - 10-^ 

Figure 3 might suggest that both correlations are produced by the same physical mech- 
anism, e.g. an increase of the temperature of the emitting gas (Trevese, Kron, & Bunone 
(2001), see also Paltani & Courvoisier (1994)). However, in order to try any comparison 
with possible models, a precise quantification of the above effects is needed. In fact, a vari- 
ety of mechanisms might in principle be responsible for a hardening of the spectrum in the 
bright phases. For instance, if AGNs are powered by supernovae explosions and variability 
is caused by a "Christmas tree" effect, then there will be an excess of blue emission in the 
bright phase (Aretxaga, Cid Fernandes & Terlevich 1997; Cid Fernandes et al. 2000). 
Even gravitational lensing (Hawkins 1996), usually thought of as achromatic, can produce a 
stronger variability in the blue than in the red band since the amplification depends on the 
size of the accretion disk which is larger in the red than in the blue band (Alexander 1995). 

Also hot spots due to instabilities of the accretion disk (Kawaguchi et al. 1998) are 
likely to enhance the blue emission. Finally the spectral energy distribution of the host 
galaxy, which is redder than the AGN and gives a stronger contribution in the faint AGN 
phase, can produce a similar effect. 



In order to quantify the spectral variations we define a spectral variability parameter 
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(SVP) representing the spectral slope changes per unit log- luminosity change: 




(4) 



where ^^(t) is the luminosity in the B band and r is an appropriate time delay. In fact we 
expect that different physical phenomena are causing different SED changes on the relevant 
time scales, which go at least from days to more than ten years. With the available light 
curves we have computed, for each QSO, l3{Tij) with Tij = ti — tj ,i,j = 1,N representing 
all the possible time differences between the N points of the light curve. Figure 4 represents 
(3{Tij) for two of the 42 QSOs considered, PG0804+762 and PG1354+213, taken as examples 
of good and poor time sampling respectively. For each bin, the mean value of the SVP is 
also shown. The uncertainty is the r.m.s. variation of the mean. No obvious trend is seen 
looking at similar plots for the entire sample: the average /3 values stay almost constant, 
within the uncertainty, at least for 1000 d. Clearly, bins at large r are less populated. An 
increase or decrease of f3 appears for some objects at larger r where, however, the statistic 
becomes poor. Apparently, an analysis of possible systematic i3 changes on time scales 3 
yr requires a time base larger than the 7 yr of the present one. For the following analysis 
we define the SVP of each QSO as the mean value (3rn in a single bin < r < 1000 d. Since 
/?(Tjj) values are not independent, the above computation of the uncertainty represents an 
overestimate of the standard deviation of P^- In Figure 6 (3^ of the 42 QSOs are reported 
versus the relevant time average a. The curves are described in the following paragraphs, 
except the dot-dashed line representing a black body, which is reported for comparison with 
the results of Trevese, Kron, & Bunonc (2001). In the case of a black body (of fixed 
area), defining x = hu/kT, with T, h, k equal to temperature, Planck and Boltzmann 
constants respectively, the spectral slope is aBsix) = 3 — xe^/{e^ — 1) and the SVP is 
{da/dT)/{d\ogBjdT) = (lnlO)[l - x/(e^ - 1)] = (3bb{x). We stress that the dotted line, 
defined by the above expression, does not represent a fit to the points (no free parameters). 
For a fixed frequency logu = \og{i'Bi'R)^ and increasing temperature, a point moves on 
the curve from left to right (Rayleigh- Jeans limit a — 2,P — 0). The "average QSO", 
(a) — —0.2 ± 1.0, {Prn) — 2.2 ± 0.9, can be represented by temperature changes of a black 
body of T 10^ K, in approximate agreement with the result of Trevese, Kron, & Bunone 
(2001), which was obtained with different QSO sample and mean rest-frame frequency. 



5. 



THE EFFECT OF THE HOST GALAXY 



An independent analysis of the same light curves, performed by Cid Fernandes et al. 
(2000) in the framework of Poissonian model of variability, implies the existence of an under- 
lying spectral component, redder than the variable one, which could be identified either with 
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the non-flaring part of the QSO spectrum or with the host galaxy considered by Romano 
& Peterson (1998). We evaluate the effect of the host galaxy through numerical simula- 
tions based on templates of the QSO and host galaxy SEDs (Trevese & Vagnetti 2001b). 
Both SEDs are derived from the atlas of normal QSO continuum spectra of Elvis et al. 
(1994). We compute a synthetic QSO+host spectrum, adding to the fixed host galaxy tem- 
plate SED the average QSO spectrum with a relative weight characterized by the parameter 
77 = log(L^/L|^), where and Lf^ are the total H band luminosities of the QSO and the 
host galaxy respectively. In the Elvis et al. (1994) sample, 77 ranges from -1 to 2. Figure 5 
shows an example of composite spectrum. We want to test (disprove) the hypothesis that the 
QSO SED maintains its shape during brightness changes, while the variation of the spectral 
shape is entirely due to the contribution of the (constant) galaxy SED. We know that the 
effect of the host galaxy should be small in the case of Giveon ct al. (1999) data since: 
i) magnitudes were computed using point spread function fitting of the images, ii) the im- 
ages were limited by an aperture with the diameter depending on the seeing conditions, but 
smaller than 3 arcsec. However we don't know the appropriate value of rj. For this reason 
we perform simulations for a range of 77 values. Variability is represented by small changes 
A77, around each 77 value, with an amplitude corresponding to a r.m.s. variability ub — 0.16 
mag in the blue band. For each synthetic spectrum, representing the QSO plus host SED 
at a given time, we compute q;(z/, t) = {d\ogLy/d\ogu)tj=tj, u = ^vbVr-, then we derive 
the SVP /9. The result, for — 3 < 77 < 3, is shown in Figure 6. To check the dependence 
of the result on QSO redshift, the same computation is repeated for the maximum redshift 
of the Giveon et al. (1999) sample, z = 0.4, and shown in the same figure. Although an 
appropriate choice of rj can reproduce the observed P oi a separately, the curves are clearly 
shifted respect to the distribution of the observational points. This means that the effect of 
the host galaxy is not sufficient to account for the observed changes of the spectral shape. 
Thus the spectral variabihty is intrinsic of the active nucleus. This also implies that the 
constant red continuum, resulting from the analysis of Cid Fernandes et al. (2000), cannot 
be identified with the host galaxy. It must be, at least in part, due to the nucleus, and can 
be identified with the spectrum of the non-flaring part of the accretion disk. 



6. CHANGES OF M AND "HOT SPOTS" 

Once the spectral variability is ascribed to the nuclear component, we can ask whether 
a change of any of the parameters defining an emission model can account for the observed 

variations of the spectral shape . We considered the accretion disk model of Siemiginowska 
et al. (1995), Fiorc ct al. (1995), corresponding to a Kerr metric and modified black body 
SED, which depend on the black hole mass M, the accretion rate M and the inclination 6 
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= ^ face-on). A grid of models has been considered for log M/M© = 7.0, 8.0, 9.0, 10.0, 
m = Mc^/Le = 0.1,0.3,0.8 (where Le is the Eddington luminosity Le = ^^^^ with the 
usual meaning of symbols) and /x = cos^ = 1, 0.75, 0.5, 0.25, 0.1. A change of M produces a 
variation of both luminosity and the SED shape. On purely theoretical grounds, we know 
that the time required for the accretion disk to reach a new equilibrium condition with a 
different M value is at least of the order of the sound crossing time t sound ~ 10'^ — 10^ d 
(Courvoisier 1991), i.e. much longer than typical variability time scales. Still, it is interesting 
to see how the spectral changes between two M states compare with the observed ones, as 
done by Tripp, Bechtold & Green (1994), and by Siemiginowska et al. (1995). In Figure 6 
the two curves on the bottom right represent (5 versus a for varying M (from right to left) 
as computed for two different values of fn and for fi = 1, from the Kerr metric, modified 
black body model of Siemiginowska et al. (1995) (their table 4). The spectral variations are 
clearly smaller, on average, than the observed ones. This means that a transition e.g. from 
a lower to a higher M regime implies a larger luminosity change for a given slope variation. 
Notice that a black body of varying T and fixed area provides larger /3 values consistent with 
the observations. Thus, (3 values computed for an increase of m would better compare with 
the case of a transition to a hotter disk with larger area. Observed slope variations of two 
objects, NGC 5548 and NGC 3783, have been compared with the predictions of an accretion 
disk model by Tripp, Bechtold & Green (1994), in the UV range. They conclude that the 
observed points lie roughly on curves of constant black hole mass, giving confidence in the 
accretion disk models. However we notice that also in this case the distribution of observed 
points is steeper than the iso-mass lines, specially in the case of NGC 3783. Our result on 
the statistical sample of 42 objects indicates that this is indeed a systematic effect. 

This result suggests that transient phenomena, like hot spots produced on the accretion 
disk by instability phenomena (Kawaguchi et al. 1998), instead of a transition to a new 
equilibrium state, may better explain the relatively large changes of the local spectral slope. 
The available models of instability phenomena do not provide a spectrum of the hot spot. 
Thus we try a simple "model" based on the addition of a black body flare to the disk 
SED, represented by the average QSO SED of Elvis et al. (1994) (shown in Figure 5). 
The free parameters are temperature Tbb, and emitting area A, while the constraints are 
the amplitude of the luminosity change (e.g. in the B band) and the relevant (3 value (or 
the relevant Aa). The solution is not univocal, given the spread of the observed j3 values. 
However, AB = 0.16 mag (corresponding to the r.m.s. variability of the sample) can be 
obtained by a hot spot of Tbb ~ 2 ■ 10^ K and A = 5 ■ 10^^ cm^, producing j3 = 3.2, or 
Tbb ~ 2 ■ 10^ K, A = 1,3 ■ 10^^ cm^, giving [3 = 2.2. This is shown by the large filled squares 
in Figure 6. A sudden heating of a fraction of the disk surface is thus capable of producing 
the observed change of the SED in the B and R bands and the relevant intra-QSO a-L 
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correlation. 



7. SUMMARY AND CONCLUSIONS 

We have performed a new analysis of the B and R hght curves of a sample of 42 PG 
QSOs made available by the Wise Observatory group (Giveon et al. 1999). We have shown 
the existence of a positive v-z correlation when the analysis is restricted to a small luminosity 
bin. 

This correlation disappears in the sample as a whole due to the interplay of v-L, L-z 
and v-z correlations and the small redhift range. The slope of the v-z relation is consistent 
with the previous findings of Di Clemente et al. (1996) and Trevese, Kron, & Bunone 
(2001), thus confirming that the v-z correlation can be entirely explained by the increase 
of variability with frequency, coupled with the increase of (rest-frame) observing frequency 
for higher redshift objects. The dependence of variability on redshift has been questioned in 
the past (see table 1 of Giveon et al. (1999) for a summary of previous variability studies). 
Wc stress that: i) it is a relatively weak effect {dv/dz ^ 0.1) so that it cannot be detected 
unless other competing effects (including v-L and z-L correlations) arc properly taken into 
account; ii) to get rid of spurious effects connected with cosmological time dilation and finite 
sampling time, variability must be quantified by an intrinsic index defined on the basis of the 
rest frame structure function. Once these prescriptions are adopted, the v-z correlations as 
measured in different optical-UV bands and different QSO samples appear consistent. This 
was already found by Di Clemente et al. (1996) and Trevese, Kron, & Bunone (2001). The 
present evidence, though marginal, is again quantitatively consistent and confirms previous 
results. 

The two-color light curves of the Giveon et al. (1999) sample allow for the first time 
the statistical analysis of SED variability of individual QSOs and the study of the spectral 
slope changes among QSOs of different luminosity, leading to the evidence of an intra-QSO 
and an inter-QSO correlation. 

We have analyzed the spectral variability by the distribution of the parameter /9(r) = 
log LB(t+T)^k)g Lflft) ^ function of the spectral slope a. We have compared with the observed 
distribution the spectral slope changes produced by the contribution of the host galaxy to 
the QSO SED, under the assumption that the nuclear spectrum maintains its shape while 
changing its brightness. We conclude that the host galaxy alone cannot be responsible for 
the observed spectral changes. Thus the spectral variation must be intrinsic of the nuclear 
component. The P-a distribution has been also compared with the SED changes of a disk 
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model for M variations. The latter appear insufficient to explain the observed spectral 
changes. Bright spots on the disk, likely produced by local instabilities, are able to represent 
the observed spectral variability. 

Since it is likely that different physical phenomena are causing variability in different 
bands and time scales, multi-frequency analyses will be ultimately needed to obtain a com- 
plete description. However even a two optical bands analysis, once performed on a statistical 
sample, provides valuable constraints on the origin of variability. This strongly suggests to 
extend the work of the Wise Observatory group both in frequency and sampling time to 
allow a more detailed comparison with possible models. 

We are grateful to the Wise Observatory Group for promptly making available their data 
to the community, and for providing us with details about the observations. We are indebted 
to Fabrizio Fiore for his help in the use of disk models and for clarifying discussions. This 
research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory, Cahfornia Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 
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Fig. 1. — Variability Sb ( 0.3±0.09 yr) versus redshift for the subsample with —23.5 < Mb < 
—22.5. The correlation coefficient r and the relevant probability (of the null hypothesis) are 
also reported. 

Fig. 2. — Variability S versus rest-frame frequency for various QSO samples (adapted from 
Di Clemente et al. (1996)). The variabihty indicator is the ensemble average of variability 
indexes over the relevant QSO samples. The frequencies of each point are defined as the 
relevant ensemble averages of the rest-frame QSO frequencies, according to the individual 
redshift and observing band. Filled symbols correspond to r = 0.3 ± 0.09 yr and open 
symbols correspond to r = 2.0 ± 0.3 yr. The stars represent the four variability indexes of 
the present analysis. 

Fig. 3. — Instantaneous spectral slope a versus monochromatic luminosity L^g . Regression 
lines a — L^^ are reported for each QSO. Dotted curves represent black bodies of different 
areas with temperature T increasing from bottom left to top right (Rayleigh- Jeans limit 
« = +2). 

Fig. 4. — The spectral variabihty parameter P{Tij) as a function of the time lag between the 
observations for PG0804-I-762 and PG1354-I-213. Filled circles represent the mean value of 
(3{Tij) in intervals of 500 d. Error bars represent the r.m.s. uncertainty of the mean. 

Fig. 5. — Synthetic spectra (thin continuous lines) resulting from the composition of the 
average QSO spectrum (dotted fines) and the host galaxy template spectrum (thick line), 
both taken from Elvis et al. (1994), rj = \og{L%/L%) = -1, 0, 1 

Fig. 6. — The spectral variability parameter f3„i (see Eq. 4) versus the average spectral slope 
for each QSO of the sample. Error bars represent the r.m.s. statistical uncertainties, which 
are mainly due to intrinsic variabihty. Dot-dashed fine: sequence of black bodies of different 
temperatures (increasing from left to right); continuous line: spectral variability due to the 
presence of the host galaxy, in the QSO rest- frame, for 77 ranging from -3 (left) to 3 (right); 
dotted line: the same as continuous line, for a QSO at the maximum redshift of the sample, 
z = 0.4; thick continuous line: spectral variability caused by an increase of m = from 
0.1 to 0.3 in the model of table 4 of Siemiginowska et al. (1995) (Kerr + modified black 
body), M increasing from 10^ to 10^'^ Mq (from right to left); thick dashed line: the same 
as above for rh increasing from 0.3 to 0.8. Large filled squares represent the hot spot model 
for T = 2 • 10^ K (upper) and T = 2 • 10^ K (lower) (see text). 
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